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ABSTRACT 

We report an observation of a small scale flux-tube that undergoes kinking 
and triggers the macrospicule and a jet on November 11, 2010 in the north 
polar corona. The small-scale flux-tube emerged well before the triggering of 
macrospicule and as the time progresses the two opposite halves of this omega 
shaped flux-tube bent transversely and approached towards each other. After ~ 
2 minutes, the two approaching halves of the kinked flux-tube touch each-other 
and internal reconnection as well as energy release takes place at the adjoining 
location and a macrospicule was launched which goes upto a height of 12 Mm. 
Plasma starts moving horizontally as well as vertically upward along with the 
onset of macrospicule and thereafter converts into a large-scale jet which goes up 
to ~ 40 Mm in the solar atmosphere with a projected speed of ~ 95 km s . We 
perform 2-D numerical simulation by considering the VAL-C initial atmospheric 
conditions to understand the physical scenario of the observed macrospicule and 
associated jet. The simulation results show that reconnection generated velocity 
pulse in the lower solar atmosphere steepens into slow shock and the cool plasma 
is driven behind it in form of macrospicule. The horizontal surface waves also 
appeared with the shock fronts at different heights, which most likely drove and 
spread the large-scale jet associated with the macrospicule. 

Subject headings: magnetohydro dynamics (MHD) : sun — chromosphere: sun - 
corona 
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1. Introduction 

Macrospicules are giant spicules, mostly observed in polar coronal holes, reaching upto 
a height between 7-45 Mm above the solar-limb with a life-time of 3-45 minutes (e.g., Bohlin 
et al. 1975; Sterling 2000; Wilhelm 2000). A number of mechanisms have been proposed for 
the formation of such plasma ejecta, e.g., gas pressure pulse (Hollweg 1982), velocity pulse 
(Suematsu et al. 1982; Murawski et al. 2011). Shibata (1982) suggested that if reconnection 
takes place in upper chromosphere/lower corona (lower chromosphere/photosphere) 
macrospicule can be triggered due to magnetic reconnection (evolution of slow shocks). 
Alternative mechanisms have also been reported for the formation of such plasma ejecta 
(e.g., Moore at al. 1977; Habbal & Gonzalez 1991; Mustsevoi & Solovev 1997; De Pontieu 
et al. 2004; Kamio et al. 2010, and references cited there) 

Apart from the spicules and macrospicules in the lower solar atmosphere, the large-scale 
jets have significant role in mass and energy transport upto the higher corona, as well as 
in destabilizing large-scale coronal magnetic fields leading the eruptions (e.g., Innes et al. 
1997; Isobe & Tripathi 2006; Chifor et al. 2006; Culhane et al. 2007; 2007; Filippov et al. 
2009; Tripathi et al. 2009; de Pontieu et al. 2010; De Pontieu et al. 2011; Judge et al. 
2012). Using STEREO/EUVI (Kaiser et al. 2008), a variety of solar jets on the basis of 
their sizes, life-times have been reported (Nistico et al. 2009). The magnetic reconnection 
was found to be one of the drivers of coronal jets (Yokoyama et al. 1995; Innes et al. 1997; 
Culhane et al. 2007; Chifor et al. 2008; Filippov et al. 2009; Nishizuka et al. 2009). Pariat 
et al. (2009) have shown that reconnection generated nonlinear Alfven waves can produce 
the polar coronal jets. Alternatively, the MHD pulse-driven models are also employed for 
the triggering of various large-scale solar jets (e.g., Srivastava & Murawski 2011; Srivastava 
et al. 2012; Kayshap et al. 2013, and references cited therein), and supported by recent 
observations (Morton et al. 2012). 
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The inter-relationship of spicules and jets is also important to understand their 
formation processes. Kamio et al. (2010) have shown the association of jet with the 
macrospicule, which was triggered by a twisted magnetic flux-rope. Close association of 
the polar surges with macrospicules is also reported by Georgakilas et al. (2001) . Moore 
et al. (2011) have suggested that granule-size emerging bi-poles (EBs) can trigger the 
longer spicules and associated Alfven waves, while the larger EBs can form X-ray jets. In 
this letter, we report firstly the evidence of the activation of a small-scale bipolar twisted 
flux-tube in the lower polar corona, which undergoes internal reconnection and triggers a 
macrospicule and associated coronal jet. We study the relationship between the formation 
of macrospicule and associated jet, as well their most likely triggering mechanism using 
SDO/AIA observations and numerical simulation. In Section 2, we discuss the observations 
of macrospicule and associated jet. In Section 3, we present the results of numerical 
simulations. We outline discussion and conclusions in last section. 



2. Observations of the Macrospicule and Associated Jet 

High resolution observations (0.6 "per pixel with cadence 12 s) from the Atmospheric 
Imaging Assembly on board the Solar Dynamic Observatory (SDO) , using its various filters 
sensitive to the plasma at different temperatures (AIA; Cheimets et al. 2009, DelZanna et 
al. 2011, Lemen et al. 2012, O'Dwyer et al. 2012), recorded the origin of macrospicule and 
associated polar jet on 11 Nov 2010 (cf., Figure 1 & MS-Jet-304.mpeg). For this study, 
we have used 304 A image sequence to study the triggering of macrospicule and jet that 
occurred at the north pole of the Sun. 

Figure 1 (top-panel) displays some selected snapshots of SDO/AIA 304 A during the 
evolution of macrospicule and associated polar jet. A bipolar, omega-shaped, small-scale 
flux-tube is emerged at the polar cap around 00:58:08 UT (top-left panel). During the 
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evolution of this flux-tube between 00:58:08-00:59:44 UT (~ 95 sec), the flux-tube shows 
some transverse bending of magnetic surfaces on its both the halves that may be the 
signature of the evolution of kink perturbations (cf., blue arrows in the snapshot 00:58:44 - 
00:58:56 UT, and cartoon in bottom panel). The combined effect of apparent rotation and 
transverse bending, generates the kinked flux-tube at smaller spatial-scales, which undergoes 
internal reconnection and further leads the macrospicule and jet (cf., MS-Jet-304.mpeg). 
To the best of our knowledge, this is the first direct observation of the evolution of kink 
perturbations in the small-scale flux-tube, which further enables internal reconnection and 
leads the formation of macrospicule and associated jet that move along ambient open field 
lines (cf., Figure 2, and cartoon in Figure 1). The two halves and legs of the flux-tube 
come closer to each other and merge to produce a brightening (cf., snapshot on 00:59:44 
UT), which is very likely a signature of internal magnetic reconnection between the two 
opposite halves of a small-scale flux-tube. A macrospicule is triggered at the same time (cf., 
01:00:56 UT snapshot in Figure 2, and cartoon in Figure 1). The macrospicule reaches up 
to a height of ~12 Mm with a speed of about 80 km s _1 , and faded within ~ 4 minute. 
During this time, a jet like feature evolved which grew in vertical as well as horizontal 
direction (cf., 01:00:56 UT - 01:02:08 UT snapshots in Figure 2). The detailed plasma 
dynamics can be seen in the associated movie MS-Jet-304.mpeg. We also investigated 
other high temperature filters of AIA (e.g. 171,211, 193 A ) where jet plasma fronts are 
clearly evident, however, the macrospicule material is only visible in 304 A channel (cf., 
MS-Jet-171.mpeg, MS-Jet-211.mpeg, MS-Jet-193.mpeg). Figure 3 displays a slit position 
along the macrospicule and jet (left panel) and corresponding height-time diagram (right 
panel). The two paths represented by pink and yellow colors respectively show the ejecting 
and down-flowing jet material. The denser core of the large-scale jet plasma goes up to 
a height ~40 Mm with a speed of ~ 95 km s _1 , while its some fainter traces are evident 
upto ~ 50-60 Mm. In the first ~10 minutes, the large-scale jet plasma material seen to be 
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moving upward and later started to fall back towards the surface within total life-time of ~ 
24 minutes. The average outflow speed and acceleration were ~ 95 km s _1 and 89.015 m 
s -2 , while the average down-flow speed and downward acceleration are respectively 154 km 
s" 1 and -649.99 m s" 2 . 

Small-scale flux-tube in the polar corona, which is most-likely imposed by the kink 
perturbations due to asymmetric bending of magnetic surfaces on its both the halves, 
creates a reconnection diffusion region on the temporal scale of ~ 95 sec. This region has 
length (2L) ~ 12 Mm, that is the height of the emerged flux-tube at 00:58:08 UT. The 
half-width (S) is ~ 3.6 Mm that is the separation between the two opposite halves of the 
kinked flux-tube (cf., snapshot on 00:58:08 UT in Figure 1). Two opposite halves of the 
small-scale flux-tube approach to each-other with a supersonic speed of ~ 76 km s _1 , (sound 
speed in choromosphere is 15 km s _1 ) and reconnected. The height of the reconnection site 
is around ~ 3.0-4.0 Mm in the polar region from the anchored footpoint of tube suggesting 
chromospheric reconnection which may trigger the evolution of slow shocks excited by the 
velocity pulse (Shibata 1982; Murawski et al. 2011; Kayshap et al. 2013) and the observed 
plasma dynamics. Considering this as basic scenario, we numerically simulate the observed 
macrospicule and jet dynamics. 



3. Numerical Model of the Macrospicule and Jet 

We perform a 2-D MHD numerical simulation to understand the physics of the observed 
macrospicule and jet using FLASH code (Lee & Deane 2009) with an assumption that the 
polar corona is gravitationally stratified. The set of equations solved using FLASH code are 
as follows: 
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dg 
di 



+ V-(f?V) = 0, 



(1) 



q-zt + q(v- V)V 



Vp + -(V xB)xB + gg, 



(2) 



<9p 



+ V • (pV) = (1 - 7 )pV ■ V 



(3) 



<9B 

a* 7 



Vx(VxB), V-B = 0. 



(4) 



Here g, V, B, p — ^gT, T, 7 = 5/3, g = (0, — g) with its value g = 274 m s~ 2 , m, fc B , 
are respectively the mass density, flow velocity, magnetic field, gas pressure, temperature, 
adiabatic index, solar gravitational acceleration, mean particle mass, and Boltzmann's 
constant. Radiative cooling as well as thermal conduction are not included in our model as 
we are only interested in the dynamics at this instance. 

3.1. Equilibrium Configuration 

We assume that the solar atmosphere is in the static equilibrium (V e = 0) with a force 
free magnetic field, 



such that it satisfies the current free condition, V x B c = 0, and it is specified by the 
magnetic flux function, A, as 




(5) 




(6) 



Here the subscript e corresponds to equilibrium quantities. 
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We set a weakly curved arcade type magnetic field configuration by choosing 

A(x, y) = B A B cos (x/A B )exp[-(y - y r )/A B ] ■ (7) 

Here, B is the magnetic field at the reference level, which is initial location of the pulse 
y = y r , and the magnetic scale- height is 

A B = 2L/tt . (8) 

We set and hold fixed L = 200 Mm and y r = 10 Mm. 

As a result of Eq. (JSJ) the pressure gradient is balanced by the gravity force, 

" Vp c + £? c g = . (9) 

With the ideal gas law and the y-component of Eq. ([9]), we arrive at 



Pe(y) = Po exp 

where 



dy 



My') 



&fe) = ^Afe)' <10) 



Hv) = k B T e (y)/(mg) (11) 

is the pressure scale-height, and po denotes the gas pressure at the reference level that we 
choose in the solar corona at y x = 10 Mm. 

We take an equilibrium temperature profile T e (z) (cf., bottom- left panel of Figure 
3) for the solar atmosphere that consists of VAL-C atmospheric model of Vernazza et al. 
(1981) and obtain the corresponding gas pressure and mass density (not shown) using 
Eq. (fit)]) . In our simulation the transition region is located at y ~ 2.7 Mm. An extended 
corona and chromosphere is considered respectively above and below this with requirement 
of temperature minimum at y ~ 0.9 Mm 
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3.2. Perturbations 



We impulsively perturb the system in equilibrium by a Gaussian velocity pulse V that 
is nearly parallel to the ambient magnetic field lines, viz., 

(x - x ) 2 + (y - y ) 2 ~ 



V\\ (x, y, t — 0) = A v exp 



w- 



(12) 



Here A v is the amplitude of the pulse, (xq, yo) is its initial position and w is its width. We 
set and hold fixed A v = 7.5 km s _1 , Xq = Mm, y = 0.9 Mm, and w = 0.2 Mm. 



3.3. Results of the Numerical Simulation 

We set the simulation box as (-10,10) Mm x (0,40) Mm and impose the boundary 
conditions by fixing all plasma quantities to their equilibrium values in time for x- and 
y-directions, while all plasma quantities remain invariant along the z-direction. In the study, 
we use AMR grid with minimum (maximum) level of refinement set to 3(8) (cf., Top-left 
panel of Figure 3). We launch the velocity pulse in the lower solar atmosphere, which is 
considered to be excited by the chromospheric reconnection generated energy between the 
opposite halves of the omega-shaped bipolar flux-tube (cf., Eq. 12). As the pulse propagates 
in gravitationally stratified atmosphere, it converts into a slow shock at higher altitudes. 
As a result, a low pressure region develops behind it and drives cool chromospheric plasma 
upwards. This lagging plasma exhibits the properties of the observed macrospicule. 

Figure 5, displays key snapshots of the simulation results as the temperature- map 
(color) and velocity (arrows) of the plasma. At 100 s (first image), the shock front of the 
initially launched pulse has reached upto ~ 6 Mm. The chromospheric plasma, however, 
lags behind the shock front and has reached up to a height of ~4 Mm. This lagging is 
primarily due to the rarefaction of the plasma behind the shock front. Rarefaction is the 
low pressure region behind the shock front and due to the low pressure region surrounding 
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plasma falls back in this region. The falling back material pushed back due to the larger 
photospheric density, which is basically lagging of plasma. By the time the chromospheric 
plasma reaches up to a height of ~10 Mm at 250 s (see the right-top panel), the shock front 
has reached a height more than 30 Mm. Our observational finding shows that macrospicule 
goes up-to ~12 Mm with its width ~4 Mm and the life-time of ~4 minutes. The simulation 
results presented here approximately mimic the observed dynamics of the macrospicule, 
which is the propagation of cool plasma. 

At t=400 s, there are two bumps present on each side of the main macrospicule, 
which are very likely created by the secondary shocks. Due to these shocks the base of the 
spicule is getting wider, i.e., the surrounding material is moving upward along with the 
macrospicule, which is most likely the spreading plasma as observed in form of the initiation 
of jet. At this moment, the central part material is suppressed compared to the bumps due 
to the dominance of the down-flowing material at the center. The bumps on each side of the 
central part goes up-to ~ 11 Mm while the central part suppressed up-to ~ 8 Mm. The next 
snapshot (at t = 450 s) , again shows that central region material is continuously falling back 
towards the solar surface suppressing it while the side wise plasma are moving upward in 
the solar atmosphere and most likely form the core of jet. In observations, the macrospicule 
fades, and the large-scale jet plasma evolves around it, which may have striking similarity 
with this model result. This quasi-periodic rise and fall with minimum time-scale of 200 s of 
the cool chromospheric material is clearly visible over the longer duration in simulation due 
to the arrival of wave trains of slow shocks (cf., Figure 4 and SIM-Jet.mpg). This time-scale 
though depends upon pulse-strength, steepening of slow shocks and their reflection from 
transition region creating down-flows. Moreover, the simulated plasma ejection creates 
horizontal surface waves at higher altitudes, which are clearly visible with comparatively 
high temperature shock fronts at different heights (cf., Figure 4 Sz Sim-Jet.mpg). These 
surface waves may interact with the various layers of the upper atmosphere upto 40 Mm 
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(cf., Fedun et al. 2011), and most likely triggers the horizontal and vertical spread of 
the multi-temperature jet (cf., MS-Jet-171.mpeg, MS-Jet-211.mpeg, MS-Jet-193.mpeg). 
Quasi-periodic rise and fall of comparatively cooler (cf., MS-Jet-304.mpeg) macrospicule 
and jet may not be evident in the observational base-line as the upcoming pulse trains may 
be observed by down-falling material (Srivastava & Murawski 2011). 

4. Discussion and Conclusions 

Using the high resolution observations of SDO/AIA at 304 A , we studied the detailed 
evolution of a macrospicule and associated jet recorded on 11 November, 2010. In addition, 
we performed numerical simulation to qualitatively match with the observed macrospicule 
and jet using VAL-C model of the solar atmosphere and FLASH code. To the best of 
our knowledge, this is the first direct evidence of formation of a macrospicule due to the 
magnetic reconnection between two opposite halves of an emerging small-scale, kinked 
bipolar loop at lower altitude in the polar corona. The observed kinked small-scale flux-tube 
may also be a rotating helical structure that undergo in internal reconnection to trigger the 
jet (e.g., Patsourakos et al. 2008; Nistico et al. 2009; Pariat et al. 2009). The triggering 
of the macrospicule was followed by the evolution of a large-scale jet. These types of 
kinked flux-rope, internal reconnection, and related plasma dynamics were only observed 
in the large-scale active regions, sometime leading to large-scale coronal eruptions, e.g. 
prominences, CMEs etc. (e.g., Torok & Kliem 2004; Tripathi et al. 2007; Tripathi et al. 
2009; Srivastava et al. 2010; Kliem et al. 2010; Srivastava et al. 2013a; Srivastava et 
al. 2013b, and references there in) . However, the analogous conditions of the evolution 
of kinked bipolar-loop at small spatio-temporal scale in the polar corona is observed for 
the first time as episodic mechanism for the triggering of macrospicule and jet. The 
macrospicule goes up-to ~ 12 Mm with the projected speed ~ 80 km s _1 , and had a life 
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time of ~4 min. The brightened and core plasma of the large-scale jet goes up-to ~ 40 
Mm in the solar atmosphere with the projected speed ~ 95 km s" 1 , and its life time was ~ 
24 minute. The standard jet models deal with the direct reconnection driven forces (jxB) 
between open and closed fields lines in the corona leading to the jet plasma propulsion 
(Yokoyama et al. 1995; Nishizuka et al. 2009; Pariat et al. 2009). However, the present new 
episodic mechanism suggests that the internal reconnection in a small-scale loop in lower 
chromosphere further generates a velocity pulse that steepens in slow-shock wave trains 
propagating through ambient open field lines and triggering the dynamics of macrospicule 
and jet. 

Depending upon the height of reconnection site inside the chromosphere and amount 
of energy release during reconnection within small-scale flux-tube (e.g., Shibata 1982; 
Murawski et al. 2011; Kayshap et al. 2013, and references therein) , it most likely generates 
the velocity pulse that further converts into a slow shock and exhibits the features of 
macrospicule and associated jet. The excitation of surface waves and the motion of the shock 
fronts and associated plasma may be responsible for the formation of the observed solar 
jet. Our numerical results, therefore, approximately and qualitatively match the observed 
plasma dynamics. We conclude that the kinking and chromospheric reconnection in the 
small-scale flux-tube can be an episodic mechanism to drive the observed macrospicule and 
associated jet via secondary consequences in form of the evolution of velocity pulse and 
associated slow shocks. 
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Fig. 1. — Top-panel: SDO/AIA 304 A (~ 10 5 K) snapshots show the evolution of a small- 
scale flux-tube in the polar region to trigger the macrospicule and associated jet. Enlarged 
view of the small-scale, omega shaped bipolar loop shows that the two opposite halves of the 
kinked flux-tube approach each other and reconnect internally (00:59:44 UT) to release the 
energy in the chromosphere. Bottom panel: A schematic showing the formation of transverse 
bending and apparent rotation in small-scale flux-tube that cause its kinking to trigger the 
internal reconnection and energy release. 
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Fig. 2. — The macrospicule triggers on 00:59:44 UT due to the reconnection between opposite 
halves of flux-tube as shown in Figure 1, and finally converts into a jet (cf., 01:02:08-01:06:08 UT 
snapshot). The comprehensive dynamics of the macrospicule and jet are shown in the attached 
movie MS-Jet-304.mpeg. The cartoon in Figure 1 also depicts the scenario of the formation of 
macrospicule and iet alone ambient open field lines. 
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Fig. 3. — Left image shows the location of the slit along the jet to measure the height-time 
profile. Right snapshot shows the height-time profile of the jet. Two different paths on 
this image show respectively the acceleration of the jet (red path) as well as its deceleration 
(yellow path). The bright and denser part of the jet goes up-to ~ 40 Mm with projected 
speed of ~ 95 km s^ 1 , while some of its fainter traces are also observed upto 50-60 Mm. 
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Fig. 4. — Top-left panel shows the grid-scheme and blocks of the simulation set-up, while 
top-right panel displays the pattern of the magnetic field vectors. Bottom-left and bottom- 
middle panels respectively demonstrate the temperature and Alfven speed profiles used in 
the numerical simulation. Bottom right most panel shows the variation of the normalized 
mass density collected at a height of 5 Mm in the transition region showing the arrival of 
shock wave trains. 
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Fig. 5. — Temperature (color maps) as well as velocity (arrows) profiles are shown at t= 100 
s, t=200 s, t=250 s, t=400 s, t=450 s, t=500. Temperature drawn in the MK while the arrows 
represent the velocity expressed in the units of 50 km s . 
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